The Quantum-Classical Transition in Nonlinear Dynamical Systems 
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Viewed as approximations to quantum mechanics, classical evolutions can violate the positive- 
semidefiniteness of the density matrix. The nature of this violation suggests a classification of 
dynamical systems based on classical-quantum correspondence; we show that this can be used to 
identify when environmental interaction (decoherence) will be unsuccessful in inducing the quantum- 
classical transition. In particular, the late-time Wigner function can become positive without any 
corresponding approach to classical dynamics. In the light of these results, we emphasize key issues 
relevant for experiments studying the quantum-classical transition. 



PACS Numbers: 03.65.Yz, 05.45.-a 

In recent years much effort has been expended, both 
theoretically and experimentally, to explore the transi- 
tion from quantum to classical behavior in a controlled 
way. In this context, the interaction of trapped, cold 
atoms with optical potentials, both time-dependent and 
independent, has become a topic of considerable interest 
and activity §]||. The experimental ability to systemat- 
ically study dissipative quantum dynamics of nonlinear 
systems is an exciting new area where the frontier be- 
tween classical and quantum mechanics may be carefully 
examined. 

In this Letter we wish to explore some of the key qual- 
itative features of the quantum-classical transition. Wc 
establish that with % fixed at a finite value, and classi- 
cal dynamical evolution equations for phase space dis- 
tribution functions viewed as approximations to the un- 
derlying quantum equations, classical Liouville and Mas- 
ter equations violate the quantum constraint of positive- 
scmidefinitcness of the density matrix: We refer to this 
property of the density matrix as 'rho-positivity'. There 
is thus a global obstruction to the classical limit aris- 
ing directly from quantum evolutions. We argue below 
that rho-positivity violation can (1) serve as a guide in 
classifying dynamical systems with regard to classical- 
quantum correspondence: weak violation as Type I, 
strong violation as Type II, and (2) explain robustness 
to decoherence in the sense of avoidance of the classi- 
cal limit as exemplified by dynamical localization in the 
(open-system) quantum delta kicked rotor (QDKR). Our 
results impact directly on the interpretation and design 
of experiments to test various aspects of the quantum- 
classical transition. 

As described in more detail below, certain nonclassi- 
cal aspects of the dynamics of the QDKR turn out to be 
stable to decohering effects of external noise and decoher- 
ence due to spontaneous emission. A dynamical descrip- 
tion in terms of the Wigner function leads to two alterna- 
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tives to explain this stability: (1) diffusion in the quan- 
tum Master equation is simply not efficient at suppressing 
quantum interference terms present in the Wigner func- 
tion, or (2) the much more intriguing possibility that the 
diffusion is successful at decohering the Wigner function, 
i.e., interference terms are suppressed and the Wigner 
function is (almost) everywhere positive, yet the late- 
time distribution is not the solution of the corresponding 
classical Fokkcr-Planck equation. We find that the sec- 
ond possibility is the one actually realized, and show how 
it arises as a consequence of the fact that the classical 
Fokker-Planck equation violates rho-positivity, while the 
quantum Master equation does not. 

The singular nature of the classical limit U — > in 
quantum mechanics has been appreciated for a long time. 
However, what has not been stressed sufficiently is the 
reason for this singular behavior, that classical dynamics 
violates unitarity and rho-positivity, and thus Ti = can- 
not be connected smoothly to % — ► 0. A simple example 
suffices to make this point clear. Let us consider as ini- 
tial condition a pure Gaussian state. Let us suppose that 
we evolve the corresponding (positive) Wigner function 
classically in some nonlinear potential (for linear systems 
classical and quantum dynamics are identical Q), then 
the distribution becomes no longer Gaussian, but is still 
positive-definite. Three possibilities now present them- 
selves: the evolved object can be interpreted as (1) a pure 
quantum state (unitarity is preserved) , (2) a mixed quan- 
tum state (rho-positivity is preserved), and (3) cannot 
be interpreted as a quantum state (rho-positivity is vio- 
lated). The first possibility can be dismissed using Hud- 
son's theorem: the only pure state with positive Wigner 
function is a Gaussian state with a (necessarily) Gaus- 
sian Wigner function But our distribution is non- 
Gaussian. As to the second, we first note that the phase 
space integral of any function of the phase space distri- 
bution is preserved under a Liouville flow. In particular 
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/ f 2 (x,p)dxdp remains constant. For Wigner functions 
this quantity is proportional to Trp 2 which is a direct 
measure of whether a state is mixed or not - since this 
measure cannot change, the evolved object is not inter- 
pretable as a mixed state. Thus we are forced to the third 
alternative, that the evolved object cannot be interpreted 
as a quantum state at all, i.e., the Weyl transform of the 
evolved classical distribution yields a 'classical density 
matrix' which is non-rho-positive j3|. 

The above analysis makes it clear that the classical Li- 
ouville equation can never arise as a formal limit of quan- 
tum theory. However, all real experiments deal with open 
systems, i.e., systems interacting with their environment, 
of which the particular case of a measuring apparatus 
(necessary to deduce classical behavior) is an important 
example. Quantum dccohcrence and conditioned evolu- 
tion arising as a consequence of such system-environment 
couplings and the act of measurement and observation 
provide a natural pathway to the classical limit as has 
been quantitatively demonstrated, e.g., in Ref. Thus, 
it is important to inquire into the role of rho-positivity 
violation in this context: When is it important, and when 
not? 

Conditions have been previously derived that apply to 
the extraction of (noisy) classical trajectories via contin- 
uous quantum measurement [||. Once these (strict) con- 
ditions are satisfied (typically in the macroscopic limit 
h <C S, where S is the system action), measurement in- 
duces classical behavior, and in this regime all systems 
are therefore Type I. However, when these conditions are 
not satisfied, as is the case in most current experiments, 
the differences are indeed important. The key point is 
that, under continuous measurement, Type I systems can 
violate the classicality conditions in the sense that indi- 
vidual classical trajectories cannot be extracted, yet ex- 
pectation values are close to the classical results, whereas 
in Type II systems, violation of the classicality condition 
also implies a violation of correspondence at the level 
of expectation values. We will demonstrate this for the 
QDKR below. 

A quantum Master equation representing a weakly 
coupled, high temperature environment often utilized in 
studies of decoherence is 
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An unraveling of the weakly-coupled, high tempera- 
ture environment, Master equation (Q) is provided by 
a continuous measurement of position. This process is 
described by a stochastic Master equation for the den- 



sity matrix p(t), conditioned on the measurement record 
(X) + £(t) with = (&nk)-V a dW/dt, § 

p(t + dt) = p-(jr[H,p]-k[X,[X,p]]) dt 

+ v / 2^k ([X,p]+-2pTrpX) dW, (4) 

where A; is a constant specifying the strength of the mea- 
surement, rj is the measurement efficiency and is a num- 
ber between and 1, and dW is a Wiener process, satis- 
fying (dW) 2 = dt. When r) = 1, the evolution preserves 
the purity of the state and can be rewritten in a way 
which allows it to be understood as a series of diffuse 
projection measurements |J. Averages over the result- 
ing Schrodinger trajectories reproduce expectation val- 
ues computed using the reduced density matrix p or the 
corresponding Wigner function fw(x,p) obtained from 
solving the Master equation (|l|). The strength of the 
measurement is related to the diffusion coefficient of Eq. 
(|) by D = h 2 k. 

When the diffusion constant D = 0, Eq. ([l]) is just the 
quantum Liouville equation for the closed system. Note 
that the linearity of the quantum Liouville equation im- 
plies that in order for the evolution to be unitary, L q 
cannot be unitary since L c i is not (the sum L r j + L q is 
unitary but not the operators separately). The familiar 
heuristic argument for obtaining the classical limit from 
the quantum Master equation is that the diffusion term 
smooths out the interference effects generated by L q in 
such a way that quantum corrections to the classical dy- 
namics are much reduced. It has also been argued, that 
at finite ft, the appropriate limiting case of the quantum 
Master equation is in fact the classical Fokkcr-Planck 
equation [setting L q = in Eqn. ([!])] rather than the 
classical Liouville equation j?j . In any case, one immedi- 
ately appreciates that if either of the classical equations 
are strongly rho-positivity-violating, i.e., are Type II then 
this implies the existence of compensatory 'large' quan- 
tum corrections in the quantum Master equation, and 
hence the above heuristic argument must fail: L q is re- 
sponsible for more than just the generation of interference 
fringes in the Wigner evolution. 

Previous work has already suggested the possibility 
that closed dynamical systems may be roughly divided 
into two types depending on the (dynamical) classical- 
quantum correspondence as follows: (1) Type I systems 
in which quantum expectation values and classical aver- 
ages track each other relatively closely as a function of 
time [[7 10 , e.g., the driven Duffing oscillator with Hamil- 
tonian, 
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and (2) Type II systems in which the quantum and clas- 
sical averages diverge sharply after some finite time, e.g., 
dynamical localization in the QDKR [jll|. The QDKR 
Hamiltonian is 
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We solved the classical and quantum Master equations 
corresponding to Eqs. (||) and (^) using a high-resolution 
spectral solver implemented on parallel supercomputers. 
The solver explicitly respects rho-positivity conservation. 

We verified that in both the QDKR and the Duffing 
oscillator numerical examples discussed below the local- 
ization condition j(J necessary to obtain classical trajec- 
tories was violated. The relevant condition is 8r/fc 3> 
(d%.F/ F) y/ d x F/2m where the force F and its derivatives 
are evaluated at typical points in phase space. For both 
cases we have in fact, 8r]k ~ (d x F/F)y / d x F/2m thus 
localization does not occur (direct numerical solution of 
the corresponding stochastic Schrodinger equation con- 
firms this result) and, as discussed earlier, a meaningful 
distinction between Types I and II is possible. As the 
value of h is reduced (with D fixed and non-zero) one 
does expect an approach to the classical limit || , though 
the trajectory in the space of D and Ti need not be simple 
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FIG. 1. Eigenvalues of the quantum density matrix (solid) 
and the classical approximation (long-dashed) computed from 
the quantum and classical Master equation evolutions for the 
QDKR at t = 6. Also shown (short-dashed) is the classical 
result for the Duffing oscillator at t = 10. 

Our numerical code returns us the classical distribution 
function, the quantum density matrix and the Wigner 
function as functions of time. We then numerically solve 
for the eigenvalues of the quantum density matrix and 
the eigenvalues of the Weyl transform of the classical 
phase space distribution (the 'classical density matrix'). 
Results of one such computation are displayed in Fig. 
1 for the QDKR (Type II) and Duffing system (Type 
I). For the QDKR, initial conditions are pure Gaussian 
Wigner functions characterized by the position width 
Ax = 2.5, momentum width Ap = 1, centered on the 
point (x,p) — (0,0), and with h = 5 and k = 10. The 
diffusion coefficient D — 0.1, corresponding to k — 0.004. 
The horizontal axis refers to the index i corresponding 
to the eigenvalues Aj, which are themselves plotted on 



the vertical axis. The solid line is a result from a nu- 
merical solution of the quantum Master equation, as ex- 
pected all eigenvalues are positive (the pure initial state 
has one eigenvalue equaling unity, the rest being zero). 
The dashed line is the corresponding result from the clas- 
sical Fokker-Planck equation, which is characterized by a 
strong contribution from negative eigenvalues. It is thus 
clear that the true quantum density matrix and that pro- 
vided by the classical approximation are in fact quite dif- 
ferent. In contrast, results from classical Duffing calcula- 
tions show a very small contribution from negative eigen- 
values [Parameter values in the particular case shown in 
Fig. 1 were m = 1, A = 10, B = 0.5, A = 10, w = 6.07, 
Ax = 0.05, Ap = 1, (x,p) = (-3, 8), h = 0.1, D = 0.02.] 
These results show how rho-positivity violation may be 
used to distinguish the two types of dynamical systems. 
An important point to emphasize is that it is sufficient 
to only carry out the classical dynamical calculation in 
order to classify a dynamical system as being Type I or 
II. (The initial condition must of course be a Wigner 
function.) Also, it should be clear that non- violation of 
rho-positivity is a necessary but not sufficient condition 
for quantum-classical correspondence in terms of agree- 
ment of expectation values. 




FIG. 2. The Wigner function negativity measure F as a 
function of time for D = and D — 0.1 for the QDKR. 

It is well-known that dynamical localization in the 
QDKR can be destroyed (in the sense that (p 2 {t)) no 
longer saturates at late times) by coupling to external 
noise or to dissipative channels (e.g., spontaneous emis- 
sion) |ll|. However, what is important to note is that 
even in the presence of quite strong coupling to these de- 
cohering channels, the evolution does not go over to the 
classical one, and in this sense the QDKR is quite dif- 
ferent from the Duffing system investigated in Ref. [Q. 
In addition to numerically solving the Master equation 
(|l|) we have investigated in detail the effects of including 
amplitude and phase noise, timing jitter in the kicked 
system, and carried out more realistic simulations taking 
into account the effects of spontaneous emission. Stabil- 
ity to decoherence in the sense above was manifest in all 



3 



of these cases. Since we have established that the QDKR 
is a Type II system (Fig. 1), this behavior is essentially 
forced: as long as the classical evolution strongly vio- 
lates rho-positivity, it is impossible for the full evolution 
to ever become close to the classical limit as the quan- 
tum corrections must always be concomitantly large. The 
question remains whether the resulting Wigner function 
at least has a classical interpretation. In order to investi- 
gate this we computed as a function of time, the quantity 
r = J dxdp(\fw \ — fw), which provides a global measure 
of negativity of the Wigner function. With D = 0, one 
sees that T increases monotonically as the Wigner func- 
tion develops the expected oscillatory structure as a con- 
sequence of quantum interference in phase space. When 
D^O, diffusion in phase space wipes out the interference 
and produces an essentially positive distribution which 
one may interpret classically. However, because rho- 
positivity must be maintained, classical evolution cannot 
connect two such positive distributions. Thus, in Type 
II systems decoherence can be successful in rendering the 
Wigner distribution positive, but yet not lead to a clas- 
sical limit. We note that in NMR systems there is an 
interesting question regarding when classical evolution of 
spins can reproduce quantum evolutions connecting spin 
states that have no entanglement and thus may be in- 
terpreted classically O). We have shown that a similar 
situation can arise even in single-particle evolution where 
entanglement is not an issue. 

Recent experiments have attempted to directly ad- 
dress the issue of environment-induced decoherence in 
the QDKR, in the context of cold atom optics jjj. De- 
spite some complications stemming from non-ideal real- 
izations, the results indicate that classical and quantum 
evolutions agree only at inordinately large noise levels. 
Note that in these experiments, parametric noise or spon- 
taneous emission was used as the decohering mechanism. 
(The non-selective Master equation for atomic motion in 
far-detuned laser light has a similar form to that of a 
particle subjected to continuous position measurement. 
However, arguments can be made that only the weak de- 
coherence regime can be accessed in this manner.) The 
parameter values in our numerical work are close to those 
actually used in the experiments thus, as with our sim- 
ulations, the experiments are not carried out in a clas- 
sical regime in the sense of Ref. ft| . And since we have 
demonstrated the strongly Type II nature of the QDKR, 
it follows immediately that to observe true classical be- 
havior, either the current experiments have to switch to 
a Type I system or employ lower values of K. Simply 
increasing the measurement constant k, or equivalcntly 
D, while it produces localization, adds noise into the tra- 
jectory which must be kept small in order to achieve the 
classical limit. This final condition requires a reduction 
in h as k is increased p|. 
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